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ABSTRACT 


We  compared  small  mammal  community  pat- 
terns between  riparian  floodplain  and  mesquite 
habitats  at  Montezuma  Castle  National  Monument, 
Yavapai  County, Arizona.  Species  diversity, abundance, 
weight  of  adult  males,  number  of  juveniles,  num- 
ber of  reproductively  active  individuals,  longevity, 
residency  status,  and  patterns  of  microhabitat  se- 
lection were  the  combined  criteria  that  we  used 
to  compare  between  communities. Although  abun- 
dances of  small  mammals  tended  to  be  higher  in  flood- 
plain  habitats,  species  diversity  was  greater  in  the 
mesquite  habitat.  Results  from  the  remaining  de- 
mographic variables  were  inconsistent  with  source- 
sink  predictions,  and  the  data  do  not  substantiate 
the  hypothesis  that  riparian  floodplain  habitats  act 
as  species  sources  whereas  mesquite  upland  areas  act 
as  sinks  or  dispersal  sites  for  small  mammals. 

Percent  cover  of  trees,  perennial  grasses,  Utter, 
and  exposed  soil  were  greater  in  the  mesquite. The 


floodplain  habitat  had  a  rockier  substrate  and  was 
structurally  more  heterogeneous  with  higher  fre- 
quency of  debris  piles.  Significant  microhabitat 
separation  among  the  three  most  abundant  small 
mammal  species  (Peromyscus  boylii,  P.  eremicus,  and 
Neotoma  albigula)  was  evident  within  both  habitat 
types.  Percent  cover  by  annual  and  perennial  grasses 
and  shrubs,  substrate,  and  frequency  of  shrubs,  trees, 
and  debris  were  significant  determinants  of  small 
mammal  distribution  within  a  habitat  type.  In 
addition,we  found  that  habitat  selectivity  by  the 
three  most  abundant  species  was  significant,  indi- 
cating that  these  species  select  a  nonrandom  subset 
of  available  habitat.  Of  the  factors  that  we  exam- 
ined, habitat  selection  and  microhabitat  separation 
are  contributing  factors  that  determine  small  mam- 
mal community  assemblages  in  riparian  floodplain 
and  mesquite  habitats  in  the  Verde  Valley  of  central 
Arizona. 
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SMALL  MAMMAL  COMMUNITIES  AT  MONTEZUMA  CASTLE  NATIONAL  MONUMENT 

CHAPTER  1 
INTRODUCTION 


Riparian  communities  of  the  Southwest  are 
distinguished  by  hydric,  mesic,  and  xeric  vegeta- 
tion associations  along  ephemeral,  intermittent,  and 
perennial  streams  (Mitsch  and  Gosselink  1993). 
Although  riparian  habitats  in  the  southwest  are  rela- 
tively limited  in  area,  they  contribute  more  biotic 
diversity  to  a  region  than  adjacent  upland  vegeta- 
tion communities  due  to  a  greater  vegetational  di- 
versity, a  perennial  water  source,  and  a  more 
complex  habitat  structure  (Szaro  1989;  Szaro  and 
Belfit  1987). 

Many  comprehensive  studies  of  avian  use  of 
riparian  habitats  have  been  conducted  in  the  west- 
ern United  States  (Finch  1986,  1988,  1989; 
Hehnke  and  Stone  1978; Knopf  1985),  while  much 
of  the  research  on  riparian  vertebrates  in  the  south- 
western deserts  has  focused  on  avian  communities 
(Carothers  et  al.  1974;  Johnson  and  Haight  1985; 
Stamp  1978;  Stauffer  and  Best  1980;  Stevens  et  al. 
1977;  Strong  and  Bock  1990;  Szaro  1980).  Hence, 
the  importance  of  riparian  habitat  to  birds,  whether 
it  is  a  hydro-,  meso-,  or  xero-riparian  habitat,  is 
well  documented.  For  instancejohnson  and  Haight 
(1985)  found  that  avian  use  of  xeroriparian  (dry, 
desert  riparian  including  washes  and  arroyos)  eco- 
systems in  the  Sonoran  Desert  was  greater  when 
compared  with  the  adjacent  upland  ecosystems 
during  all  seasons  Avian  species  diversity  and  popu- 
lation densities  were  approximately  five  to  ten  times 
those  of  the  surrounding  uplands. 

The  importance  of  southwestern  riparian  habi- 
tat to  small  mammal  communities,  however,  is  not 
as  well  documented.  Although  a  plethora  of  scien- 
tific literature  exists  on  small  mammals  in  a  variety 
of  desert  ecosystems  and  on  the  mechanisms  de- 
termining their  community  composition  (Bow- 
ers 1982,  1988;  Bowers  and  Brown  1982;  Brown 
1989;M'Closkey  1978, 1981;  Morton  et  al.  1994; 
Munger  et  al.  1983;  Murray  et  al.  1986;  Price  and 
Brown  1983;Thompson  1982),  most  of  these  stud- 
ies have  focused  on  granivorous  rodents  because 
the  relatively  stable,  xeric  habitats  provide  excel- 
lent systems  for  assessing  the  effects  of  food  limita- 
tion, microhabitat  utilization,  and  interspecific 
competition  on  a  group  of  ecologically  similar  spe- 


cies (Brown  and  Munger  1985).  Of  existing  stud- 
ies focusing  on  small  mammals  in  riparian  habitat, 
most  researchers  have  compared  riparian  small 
mammal  communities  to  adjacent  upland  small 
mammal  communities,  and  found  an  increase  in 
species  diversity,  abundance,  and  reproductive  ac- 
tivity in  riparian  versus  upland  areas  (Cross  1985; 
Doyle  1990;  Geier  and  Best  1980;  Olson  and  Knopf 
1988;  Stamp  and  Ohmart  1979).  Conversely, Boeer 
and  Schmidly  (1977)  found  that  riparian  commu- 
nities in  Big  Bend  National  Park, Texas,  had  lower 
species  richness  and  diversity  indices  of  small  mam- 
mals than  the  other  major  plant  communities. 
Hence,  no  consistent  generalization  exists  regard- 
ing small  mammal  assemblages  and  diversity  in  ri- 
parian versus  adjacent  habitats. 

The  purpose  of  this  research  effort  was  to  ex- 
amine small  mammal  community  patterns  between 
riparian  floodplain  habitat  and  the  adjacent,  tran- 
sitional, mesquite  habitat.  Our  first  objective  was 
to  test  the  hypothesis  that  riparian  floodplain  at 
Montezuma  Castle  National  Monument  is  a  source 
or  is  superior  habitat  for  most  of  the  small  mam- 
mal species,  whereas  the  adjacent  mesquite  habitat 
is  a  dispersal  sink  for  these  species.  To  investigate 
this  hypothesis,  we  evaluated  the  differences  be- 
tween these  two  habitats  using  the  following  cri- 
teria: species  diversity  and  abundance,  weight  of 
adult  males,  number  of  captured  juveniles  and 
subadults,  number  of  reproductive  adults  captured, 
then  longevity,  and  residency.  We  also  examined 
dispersal  between  habitats. 

We  made  several  predictions  about  how  each 
of  the  above  criteria  would  differ  between  the  two 
habitats.  Assuming  the  floodplain  is  a  source  habi- 
tat, we  predicted  that  diversity  and  abundance 
would  be  higher  in  the  floodplain,  individual  adult 
males  would  weigh  more  in  the  floodplain,  and 
the  number  of  postbreeding  season  juveniles  would 
increase  in  the  mesquite  due  to  competitive  ex- 
clusion from  the  more  optimal  floodplain  habitat. 
Number  of  reproductive  individuals  would  be  sig- 
nificantly greater  in  the  floodplain,  and  survival  of 
individuals  would  be  greatest  in  the  floodplain.We 
also  expected  mesquite  habitat  to  have  a  higher 
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proportion  of  transient  animals.  Finally,  significant  mammals  do  not  differ  between  riparian  flood- 
dispersal  from  floodplain  habitat  (emigration)  into  plain  and  mesquite  habitat  types.  We  specifically 
mesquite  habitat  (immigration)  should  occur.  examined  whether  microhabitat  separation  and 
Our  second  objective  involved  testing  the  hy-  selection  were  both  occurring  among  small  mam- 
pothesis  that  patterns  of  microhabitat  use  by  small  mal  species  in  the  two  habitats. 


SMALL  MAMMAL  COMMUNITIES  AT  MONTEZUMA  CASTLE  NATIONAL  MONUMENT 

CHAPTER  2 
METHODS 


Study  Site 

The  study  site  was  •within  Montezuma  Castle 
National  Monument  located  in  the  Verde  Valley, 
Yavapai  County,  approximately  48.3  km  (30  miles) 
from  Cottonwood,  in  central  Arizona.This  Monu- 
ment is  approximately  344  hectares  in  size,  and 
located  in  a  desert  scrubland  transitional  life  zone 
at  the  base  of  the  Mogollon  Rim,  the  escarpment 
representing  the  southern  limit  of  the  Colorado 
Plateau. 

Vegetation  of  the  riparian  floodplain  within  the 
Monument  is  dominated  by  Chilopsis  linearis  (desert 
willow),  Fraxinus pennsylvanica  (velvet  ash),  Platanus 
wrightii  (Arizona  sycamore),  Prosopis  velutina  (vel- 
vet mesquite),  and  interspersed  characteristic 
"upland"  species  such  as  Acacia  greggii  (catclaw 
acacia),  Berberis  haematocarpa  (red  barberry  or 


mahonia),  and  Juniperus  monosperma  (one-seeded 
juniper).  This  active  floodplain  has  a  rocky  sub- 
strate, a  well-developed  channel  morphology,  and 
a  regular  flood  regime  (Figure  1). 

The  adjacent  mesquite  habitat  is  located  on  a 
higher  riparian  terrace  distinct  from  the  floodplain 
habitat.  The  vegetation  of  the  mesquite  zone  is 
dominated  by  Prosopis  velutina,  Acacia  greggii, 
Gutierrezia  sarothrae  (broom  snakeweed),  Atriplex 
canescens  (four- wing  saltbush),  Berberis  haematocarpa, 
and  several  species  of  annual  and  perennial 
grasses.  The  upland  substrate  is  a  sandy-loam  soil 
interspersed  with  rock  and  gravel  patches.  Further 
distance  from  the  creek  and  at  higher  elevations, 
Larrea  divaricata  (creosote)  associations  and  Canotia 
holocantha  (crucifixion-thorn)  associations 
develop. 


700 


Oct    Oct     Oct 
1972  1973    1974 
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Figure  1 .  Time  series  of  average  Cubic  Feet  per  Second  (CFS)  discharge  data  for  Wet  Beaver  Creek 
from  October  1972  throught  September  1993.  Data  were  not  collected  from  the  gauge  during  October 
1982  through  September  1991.  Data  were  collected  by  United  States  Geological  Survey  steam  gauge/ 
weather  station  above  Montezuma  Casde  National  Monument  near  Rimrock,  Arizona. 


TECHNICAL  REPORT  NPS/NAUMOCA/NRTR-96/11 


Live  Trapping  Methodologies 

We  chose  two  locations  within  Montezuma 
Castle  National  Monument  to  establish  live-trapping 
grids,  randomly  selecting  two  locations  from  avail- 
able riparian  floodplain  strips  with  adjacent  mes- 
quite  habitat  using  aerial  photographs.  The 
floodplain  grids  were  situated  adjacent  to  Wet  Bea- 
ver Creek  about  20-25  m  from  the  stream  bank 
with  the  rows  running  parallel  and  columns  run- 
ning at  right  angles  to  the  bank.  We  placed  mesquite 
grids  50  m  above  and  parallel  to  the  floodplain 
grids  (Figure  2).  Grids  contained  100  trap  stations 
spaced  10  m  apart  (90  by  90  meters  or  0.81  hectares). 
Each  grid  encompassed  a  distinct  habitat  type.  We 
employed  this  particular  study  design  to  determine 
if  significant  animal  dispersal  from  one  habitat  to 
another  was  occurring  between  sampling  periods. 


We  trapped  both  a  floodplain  and  a  mesquite 
grid  simultaneously  at  each  location  in  order  to 
uniformly  distribute  environmental  variables  across 
habitats.  This  design  would  also  potentially  mini- 
mize the  confounding  factors  associated  with  time. 
To  accomplish  this,  we  utilized  two  trapping 
configurations:  1)  for  the  first  half  of  the  session, 
we  set  50  traps  in  every  other  column;  2)  then  we 
moved  these  traps  to  alternate  columns  during  the 
second  half  of  the  session.  We  set  one  Sherman 
live  trap  (8  by  9  by  23  cm)  baited  with  whole  oats 
at  each  station  before  sunset  and  checked  traps  at 
sunrise  the  following  morning,  sampling  each 
location  seasonally  from  4  to  8  nights,  or  until 
adequate  recaptures  were  attained.  Up  to  8  night 
sampling  sessions  were  chosen  following  White 
et  al.'s  (1982)  guidelines. 
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Figure  2.  Montezuma  Casde  National  Monument  and  the  location  of  the  two  study  locations  in  rela- 
tion to  Wet  Beaver  Creek.  Each  location  contains  a  riparian  floodplain  and  a  mesquite  live-trapping  grid 
0.81  hectares  in  size. 
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Information  that  we  recorded  at  the  time  of 
capture  included  date,  grid  number,  trap  station, 
species,  weight  in  grams,  and  an  individual  mark 
or  number.  We  also  assessed  sex,  age,  and  repro- 
ductive condition  at  time  of  capture.  Peromyscus 
spp.,  Perognathus  inter medius,  and  Reithrodontomys 
megalotis  were  marked  with  toe-clips,  whereas  the 
larger  Neotoma  spp.  and  Dipodomys  ordii  were  marked 
using  Monel  ear-tags  (Korn  1987;  Kumar  1979). 
All  field  techniques  followed  the  American  Soci- 
ety of  Mammalogists  (1987)  recommendations  for 
acceptable  field  methods  in  mammalogical  research. 

Species  Diversity 
We  calculated  species  diversity  for  both  flood- 
plain  and  mesquite  habitats  by  jackknifmg  the 
Shannon  diversity  index  (Krebs  1989;  Magurran 
1988).  Numbers  of  individuals  of  each  species  cap- 
tured for  each  of  the  eight  trapping  occasions  were 
used  in  the  jackknife  procedure  to  produce  eight 
pseudo-values.  We  constructed  95%  confidence 
intervals  using  the  average  of  these  eight  pseudo- 
values  (Adams  and  McCune  1979). 

Abundances 

We  used  the  program  CAPTURE  to  screen 
for  uniform  trap  responses  and/or  heterogeneity 
in  capture  probabilities  for  the  three  most  abun- 
dant species:  Peromyscus  boylii  (brush  mouse),  P. 
eremicus  (cactus  mouse),  and  Neotoma  albigula 
(white-throated  woodrat)  (White  et  al.  1982). The 
program  CAPTURE  also  calculates  the  mean 
maximum  distance  moved  (MMDM)  by  animals 
captured  two  or  more  times.  We  used  this  infor- 
mation to  estimate  the  effective  trapping  area  for 
the  three  species  by  adding  a  boundary  strip  to  the 
actual  grid  area.  The  boundary  strip  was  equal  to 
the  MMDM. This  boundary  strip  compensates  for 
"edge  effect"  in  estimating  the  area  of  a  trapping 
grid  (Wilson  and  Anderson  1985). 

We  used  Chapman's  unbiased  version  of  the 
Lincoln-Petersen  estimator  to  determine  small 
mammal  population  size  and  its  associated  vari- 
ance (Krebs  1989;  Menkens  and  Anderson  1988). 
Since  the  Lincoln-Petersen  estimator  uses  data  from 
only  two  trap  periods,  data  matrices  were  con- 


densed from  four  and  eight  trapping  periods  to 
two  periods. The  first  half  of  the  trapnights  during 
a  session  was  denoted  n;  the  second  half  of  the 
trapnights,  n  .  The  number  of  unique  animals 
trapped  in  each  of  these  periods,  n  and  n  ,  and  the 
number  trapped  in  the  second  period  previously 
caught  and  marked  in  the  first  period,  m  ,  were 
then  computed.  We  constructed  95%  confidence 
intervals  around  population  estimates  for  each  trap- 
ping session  using  the  Poisson  approximation  to 
the  Binomial  Interval,  since  the  total  number  of 
individuals  captured  in  the  second  period  was  less 
than  50  (Krebs  1989). 

Five  indices  were  employed  to  evaluate  habitat 
quality:  1)  body  weight  of  adult  males;  2)  age  class; 
3)  reproductive  condition;  4)  longevity;  and  5)  resi- 
dency. 

Body  Weights 

We  averaged  body  weights  of  adult  males  within 
a  habitat  type  and  compared  the  averages  using 
independent  f-tests  (Zar  1984:126-131).  Adult  fe- 
male body  weights  were  excluded  to  eliminate  bias 
from  undetected  pregnancies. 

Age  Classes 
The  age  class  of  each  animal  was  recorded  into 
one  of  three  categories:  adult,  subadult,  and  juve- 
nile. The  criteria  that  we  used  for  aging  included 
body  size  and  color  of  pelage.  Number  of  indi- 
viduals within  each  age  class  and  within  each  habitat 
type  were  determined  and  proportions  compared 
within  a  breeding  season  (April  -  July)  using  the 
G-test,  or  the  log  likelihood  ratio  chi-square  test 
for  homogeneity  of  proportions  (Sokal  and  Rohlf 
1981:695-698). 

Reproductive  Conditions 

Reproductive  condition  of  each  animal  was 
determined  as  follows:  females  were  considered 
reproductive  if  pregnant,  lactating,  or  if  they  had 
enlarged  nipples;  males  were  considered  reproduc- 
tive if  they  exhibited  scrotal  testes  or  an  obvious 
epididymal  bulge.  We  calculated  number  of  indi- 
viduals reproductively  active  during  the  breeding 
season  (April  -  July,  both  years)  for  each  habitat 
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and  each  species  and  compared  between  habitats 
using  G-tests. 

Longevity  and  Residency 

We  expected  animals  to  survive  longer  in  the 
source  habitat.  To  investigate  if  animals  survived 
longer  in  floodplain  than  in  mesquite,  and  whether 
there  was  a  higher  proportion  of  transients  in  the 
mesquite,  we  considered  animals  as  transients  if 
captured  only  once  within  a  habitat.  As  an  indica- 
tion of  survival  or  longevity  differences  between 
habitat  types,  we  compared  proportions  of  indi- 
viduals captured  in  more  than  one  trapping  occa- 
sion with  proportion  of  transients,  using  G-tests 
(Lebreton  et  al.  1992).  Proportion  of  residents  (in- 
dividuals recaptured  within  a  trapping  occasion  and 
between  occasions)  were  also  compared  to  the  pro- 
portion of  transients  using  G-tests. 

Habitat 

We  measured  habitats  during  July  and  August 
1994,  a  period  in  which  foliage  and  cover  was 
maximal,  measuring  habitat  characteristics  within 
a  2.5-meter  radius  plot  centered  at  each  trapping 
station  on  each  grid  (n=400).We  used  a  5-meter 
line  intercept  (diameter  of  the  plot)  as  suggested 
by  Rowlands  (1994)  and  Mueller-Dombois  and 
Ellenberg  (1974).  Each  5-meter  line  intercept  was 
randomly  oriented  and  centered  on  the  Sherman 
live-trap  position.  We  recorded  length  in  cm  of 
the  line  intercepted  by  each  individual  plant  por- 
tion; we  also  recorded  rocks,  gravel,  sand,  bare  soil, 
and  Utter.  We  calculated  mean  percent  cover  by 
species  of  plant  and  mean  percent  substrate  for  each 
grid.  In  order  to  simplify  interpretation  by  reduc- 
ing the  data  set  and  eliminating  intercorrelation 
among  habitat  variables,  we  combined  the  percent 
cover  values  of  species  with  the  same  morphology 
to  form  five  vegetative  strata  used  for  all  subse- 
quent data  analyses  (annual  grasses,  perennial 
grasses,  forbs,  shrubs,  and  trees). 

We  also  recorded  physical  habitat  characteris- 
tics deemed  important  to  small  mammals  within 
the  2.5-m  radius  circular  plot;  these  included  num- 
ber of  woody  debris  piles  (debris  piles  are  defined 
as  woody  ground  structures  exceeding  10  cm  in 


depth  and  width  and  distinct  from  the  surround- 
ing ground  litter),  as  well  as  the  number  and  spe- 
cies of  trees  and  shrubs.  Dead  wood  smaller  than  3 
cm  in  diameter  was  classified  as  litter,  while  larger 
pieces  were  termed  debris. 

Vegetation  and  structural  characteristics  (veg- 
etative cover  by  strata,  substrate  cover,  and  num- 
ber of  woody  debris  piles)  were  compared  between 
floodplain  and  mesquite  using  the  Wilcoxon  two- 
sample  test  (Conover  1980:280-283).  Species  rich- 
ness or  the  total  number  of  plant  species  was 
calculated  for  each  habitat  (Magurran  1988). 

We  used  canonical  correspondence  analysis 
(CCA)  to  compare  patterns  of  microhabitat  distri- 
bution among  the  three  most  abundant  small  mam- 
mal species  (P  boylii,  P.  eremicus,  N.  albigula).  The 
program  CANOCO  (Canonical  Community  Or- 
dination) was  used  to  supply: 

1)  scores  of  species  of  small  mammals  on  the 
ordination  axis; 

2)  biplot  scores  of  the  environmental  variables 
and  centroids  for  nominal  environmental  vari- 
ables for  the  ordination  diagram; 

3)  eigenvalues; 

4)  species-environment  correlations;  and, 

5)  cumulative  percentage  variance  for  species 
data  and  for  species-environment  relation  (Ter 
Braak  1986). 

To  investigate  whether  the  observed  differences  in 
microhabitat  utilization  by  the  three  species  shown 
in  the  resulting  ordination  diagram  could  be  ac- 
counted for  by  pure  chance,  we  used  a  Monte  Carlo 
permutation  test  of  the  trap  stations  with  the  first 
eigenvalue  as  the  test  statistic.  This  permutation 
exercise  compares  randomly  calculated  eigenval- 
ues to  the  eigenvalue  calculated  from  the  data  and 
tests  for  significance  of  the  first  canonical  ordina- 
tion axis  to  the  species  distributions  (Ter  Braak 
1988). 

We  utilized  a  one-way  multivariate  analysis  of 
variance  (1-way  MANOVAJohnson  andWichern 
1992:246-248)  to  distinguish  between  microhabi- 
tats  that  were  used  by  a  species  and  microhabitats 
where  animals  were  not  captured.  For  each 
species,we  compared  habitat  variables  between  traps 
where  animals  were  captured  (used  habitat)  to  a 
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random  sample  of  traps  where  animals  were  not 
captured  (available  habitat)  (Block  and  Brennan 
1993).  Habitat  variables  included  the  eight  struc- 
tural characteristics  (frequency  of  debris  piles,  trees, 
shrubs,  and  five  substrate  variables)  and  percentage 
cover  by  the  five  vegetative  strata  also  used  in  the 
canonical  correspondence  analysis.  For  the  1-way 
MANOVA,  presence  or  absence  of  an  animal  were 
considered  the  levels  encountered,  and  the  13  habi- 
tat variables  were  independent  variables.  We  then 


ran  univariate  F-tests  to  determine  which  habitat 
variables  were  significantly  difFerent  between  used 
and  available  habitat  (Snedecor  and  Cochran  1989: 
223-224).  Frequency  of  debris  piles,  shrubs,  and 
trees  were  also  compared  for  each  species  between 
used  and  available  habitat  with  G-tests.The  above 
analyses  were  conducted  both  between  habitats 
(400  trap  stations)  and  within  a  habitat  (200  trap 
stations).  We  used  a  significance  level  of  0.05 
throughout  all  analyses. 
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CHAPTER  3 
RESULTS 


Live  Trapping  Results 

A  total  of  4290  trap-nights,  2145  in  riparian 
floodplain  and  2145  in  upland  mesquite,  were  com- 
pleted during  the  course  of  this  study.  After  ad- 
justing trap-nights  for  traps  found  closed,  but  empty, 
we  captured  a  total  of  231  individuals  520  times. 
Recapture  rates  averaged  72.0%,  and  for  most  trap- 
ping occasions,  recapture  rates  exceeded  75%,  sug- 
gesting that  each  grid  location  was  sampled 
adequately  (Table  1). 


We  captured  1 0  species  representing  three  fami- 
lies of  rodents  (Table  2).  A  total  of  seven  species 
were  captured  in  the  floodplains,  and  nine  species 
in  the  mesquite  habitats.  We  captured  three 
cricetine  species  (Muridae:  Cricetinae)  most  fre- 
quently in  both  floodplain  and  mesquite  habitats: 
P.  boylii  (66  first-captures,  181  total  captures),  P. 
eremicus  (65  first-captures,  145  total  captures),  and 
N.  albigula  (55  first-captures,  118  total  captures). 
Because  other  species  generally  had  insufficient 


Table  1.  Summary  of  live-trapping  occasions,  dates,  adjusted  effort  (trap-nights  adjusted  for  traps  closed  before 
entry), and  recapture  rate  (calculated  using  all  species  combined)  at  Montezuma  Casde  National  Monument  from 
April  1993  -  May  1994. 


Adjusted 

Grid 

Dates 

Trap-nights 

Recapture 

Location 

(effort) 

rate  (total) 

Floodplain  I 

April  9-11,  1993 

100 

50.0% 

April  23-25 

1993 

99 

100.0% 

Mesquite  I 

April  9-11 

1993 

100 

100.0% 

April  23-25 

1993 

93 

80.0% 

Floodplain  II 

June  1-9 

1993 

399 

52.4% 

Mesquite  II 

June  1-9 

1993 

399 

64.7% 

Floodplain  I 

July  6-14 

1993 

392 

88.6% 

Mesquite  I 

July  6-14 

1993 

385 

80.0% 

Floodplain  II 

January  8-12 

1994 

192 

54.5% 

Mesquite  II 

January  8-12 

1994 

195 

66.7% 

Floodplain  I 

April  8-10 

1994 

94 

75.0% 

April  12-14 

1994 

96 

72.0% 

Mesquite  I 

April  8-10 

1994 

95 

75.0% 

April  12-14 

1994 

93 

55.6% 

Floodplain  II 

May  17-25 

1994 

387 

80.0% 

Mesquite  II 

May  17-25 

1994 

396 

50.0% 

Floodplain  I 

June  17-25 

1994 

386 

80.0% 

Mesquite  I 

June  17-25 

1994 

389 

71.4% 

Total  adjusted  trap-nights 

4290 

Average  recapture  rate 

72.0% 
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Table  2.  Numbers  of  small  mammal  captures  by  species  and  habitat  type  at  Montezuma  Castle  National 
Monument  from  April  1993  through  June  1994.  Numbers  in  parentheses  are  first  captures  within  a 
trapping  occasion. 


Species 


Number  of  Individuals 


Floodplain 


Mesquite 


Thomomys  bottae 
(Botta  s  pocket  gopher) 

Dipodomys  ordii 
(Ord's  kangaroo  rat) 

Perognathus  intermedius 
(Rock  pocket  mouse) 


0(0) 
0(0) 
2(2) 


2(2) 
15(7) 
0(0) 


Reithrodontomys  megalotis 
(Western  harvest  mouse) 

Peromyscus  eremicus 
(Cactus  mouse) 


1(1) 
76  (34) 


8(6) 
69  (31) 


Peromyscus  maniculatus 
(Deer  mouse) 


0(0) 


9(3) 


Peromyscus  spp.  (juveniles) 


Peromyscus  boylii 
(Brush  mouse) 


17(9) 
159  (54) 


9(8) 
22  (12) 


Neotoma  albigula 
(White-throated  woodrat) 

Neotoma  stephensi 
(Stephen's  woodrat) 

Neotoma  mexicana 
(Mexican  woodrat) 

Totals 


71  (31) 

0(0) 

11(5) 

337  (136) 


47  (24) 
2(2) 

(0) 
183  (95) 
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captures  for  quantitative  comparisons,  except  for 
calculations  of  species  diversity,  we  restricted  the 
remaining  analyses  to  these  three  species.  Effective 
trapping  area  for  P.  boylii  was  1.06  hectares,  0.99 
for  P.  eremicus,  and  1.02  for  N.  albigula. 

Species  Diversity 

Species  diversity  was  significantly  greatest  in 
the  mesquite  habitat  (Figure  3).  Shannon  diversity 
indices  averaged  1 .30  (+0.05  in  standard  error)  for 
floodplain  habitat  and  1.78  (+0.15  in  standard  er- 
ror) for  mesquite.  This  did  not  follow  our  initial 
prediction  of  greater  diversity  in  the  floodplain. 

Abundances 

For  most  species  and  most  trapping  occasions, 
CAPTURE  chose  the  model  M  ,  constant  cap- 
ture probabilities  and  no  evidence  of  heterogeneous 
trap  responses,  or  trap  shy  or  trap  happy  animals 
(Appendix  A).  Capture  probabilities  (p-hats)  were 


calculated  by  CAPTURE  depending  on  the  model 
chosen.  If  the  model  was  M  ,  or  there  was  evi- 
dence of  a  heterogeneous  trap  response,  then  p- 
hat  was  calculated  by  averaging  capture 
probabilities.  In  order  to  adequately  estimate  popu- 
lation size  for  an  area  when  N<100  using  the  pro- 
gram CAPTURE,  the  average  probability  of 
capture  (p)  must  be  in  the  neighborhood  of  0.5. 
Since  most  p-hats  were  less  than  0.30,  the  Lin- 
coln-Petersen  estimator  was  used  to  calculate  popu- 
lation size  for  all  trapping  occasions. 

Abundance  of  small  mammals  differed  between 
floodplain  and  mesquite  habitats  (Figure  4). Abun- 
dance for  each  species  was  generally  greater  in 
floodplain,  but  not  significantly  so  except  for  P. 
boylii.  Peromyscus  eremicus  abundance  was  generally 
greater  in  floodplain  except  during  June  1993 
where  6  individuals  were  captured  in  the  mesquite, 
and  none  in  the  floodplain.  Neotoma  albigula  were 
also  more  abundant  in  the  floodplain  .Values  miss- 
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Figure  3.  Comparison  of  small  mammal  diversities  in  riparian  floodplain  and  mesquite  habitats  at 
Montezuma  Casde  National  Monument  from  April  1993  -June  1994.  Bar  values  are  the  averages  of  all 
pseudovalues  from  jackknifing  the  Shannon  diversity  index.  Lines  on  bars  are  95%  confidence  intervals. 
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Figure  4.  Lincoln-Petersen  population  estimates  for  three  species  of  small  mammals  {Peromyscus  boylii,  P.  eremicus, 
and  Neotoma  albigula)  at  Montezume  Castle  National  Monument,  Yavapai  County,  Arizona  from  April  1993  - 
June  1994.  Bars  are  95%  confidence  intervals. 
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ing  in  Figure  4  indicate  that  sample  sizes  were  in- 
sufficient to  calculate  reliable  Lincoln-Petersen  es- 
timates. 

Body  weights 

We  found  that  average  weight  of  adult  males 
differed  significantly  between  habitat  types  only  for 
N.  albigula  (Figure  5).  Larger-sized  males  of  this 
species  were  found  in  the  mesquite  which,  again, 
did  not  follow  our  initial  prediction. 

Age  Classes 

In  the  nonbreeding  season  (January  1994),  we 
captured  only  six  subadults  (five  P.  eremicus  and  one 
N.  albigula).  No  juvenile  animals  were  captured. 
These  sample  sizes  were  insufficient  to  analyze  age 
proportions  for  the  nonbreeding  season.  During 
the  breeding  season  (April  -June, both  years), ju- 
venile and  subadult  age  classes  were  combined  to 
increase  sample  sizes  within  cells  for  the  G-test. 
Proportion  of  juveniles  and  subadults  during  the 
breeding  season  for  each  species  did  not  differ  sig- 
nificantly between  the  two  habitats  except  for  N. 
albigula  (Figure  6).  Eight  subadults  of  this  species 
were  captured  in  the  floodplain  habitat  and  none 
in  the  mesquite  (G=  10.46,  d.f  =  l,p<0.05). 

Reproductive  Conditions 

We  captured  no  animals  in  reproductive  con- 
dition during  winter  trapping  (January  1994). 
During  the  breeding  season,  a  significantly  larger 
proportion  of  P.  eremicus  individuals  were  repro- 
ductively  active  in  the  mesquite  in  proportion  to 
nonreproductive  individuals  (Figure  7). There  was 
no  significant  difference  in  breeding  activity  be- 
tween floodplain  and  mesquite  for  the  other  two 
species. 

Longevity  and  Residency 

We  expected  individuals  of  a  species  to  live 
longer  in  the  floodplain,  if  the  floodplain  habitat  is 
a  source  habitat  for  a  species.  Proportion  of  ani- 
mals surviving  between  trapping  occasions  did  not 
differ  significantly  between  habitat  types  for  any 
species,  though.  We  also  expected  the  mesquite 
habitat  would  show  a  higher  proportion  of  tran- 


sient animals  (those  animals  captured  only  once). 
This  was  true  only  for  P.  boylii  (Figure  8)  .We  found 
significantly  more  transients  in  proportion  to  resi- 
dents in  the  mesquite  for  this  species  (G=4.837, 
d.f.=l,p<0.05). 

Dispersal 

We  observed  few  instances  of  large-scale  dis- 
persal of  animals  between  habitats  over  the  two 
years  of  this  study.  Three  animals,  all  P.  boylii,  moved 
between  floodplain  and  mesquite  grids  in  1994; 
one  adult  male  moved  from  mesquite  into  flood- 
plain  at  location  1,  one  adult  female  moved  from 
floodplain  into  mesquite  at  location  1,  and  one 
adult  female  moved  from  mesquite  to  floodplain 
at  location  2. 

Habitat  Results 

We  found  significantly  greater  percent  cover 
of  trees,  perennial  grasses,  and  litter  in  mesquite 
than  in  the  floodplain  (Table  3).  Debris  pile  fre- 
quency and  percent  cover  of  herbaceous  vegeta- 
tion were  higher  for  floodplain  habitats  than 
mesquite.  Because  floodplain  grids  included  por- 
tions of  stream  channels,  we  found  a  greater  per- 
cent rock  substrate  in  floodplain  habitats  whereas 
we  found  percent  exposed  soil  was  higher  in  the 
mesquite. 

We  also  found  a  higher  plant  species  richness 
in  floodplain  habitats  versus  the  mesquite  (Table 
4).  A  total  of  29  species  were  encountered  in  the 
floodplain  and  only  21  species  in  the  mesquite. 
For  the  floodplain,  the  most  frequent  shrubs  were 
Chilopsis  linearis,  Gutierrezia  sarothrae,  Acacia  greggii, 
Berberis  haematocarpa,  and  Brickellia  californica 
(pachaba). Predominant  trees  were  Platanus  wrightii, 
Juniperus  monosperma,  Prosopis  velutina,  Fraxinus 
pennsylvanica,  and  Juglans  major  (Arizona  walnut). 
A  total  of  four  species  of  perennial  grasses,  several 
species  of  annual  grasses,  six  species  of  herbs,  and 
nine  different  shrub  species  were  found  in  flood- 
plain  habitat. 

In  mesquite  habitat,  predominant  shrubs  were 
A.  greggii,  G.  sarothrae,  Atriplex  canescens,  Marrubium 
vulgare  (common  horehound),and  B.  haematocarpa. 
Tree  species  were  J.  monosperma  and  P.  velutina.  A 
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Figure  5.  Comparison  of  average  weight  in  grams  of  adult  male  P.  boylii,  P.  eremicus,  and  N.  albigula  (both  years 
combined)  between  floodplain  and  mesquite  habitatat  Montezuma  Castle  National  Monument,Yavapai  County, 
Arizona.  Numbers  above  bars  are  sample  sizes  (first  captures).  Lines  on  bars  are  standard  errors. 
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Figure  6.  Breeding  season  age  compostition  (proportion  of  adults  compared  to  juveniles  and  subadults)  between 
floodplain  and  mesquite  habitats  at  Montezuma  Castle  National  Monument,  Yavapai  County,  Arizona  during 
April  1993  -June  1994  excludingjanuary  1994  trapping  occasion.  *Log  likelihood  ratio  chi-square  value=10.464, 
d.f.=l,/K0.05. 
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Figure  7.  Proportion  of  adults  in  reproductive  condition  during  the  breeding  season  for  three  species  of  small 
mammals  at  Montezuma  Casde  National  Monument,  Yavapai  County,  Arizona  during  April  1993  -June  1994 
excluding  January  1994  trapping  occasion.  *Log  likelihood  ratio  chi-quare  value=5.104,  d.f.=l,/K0.05. 
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Figure  8.  Proportion  of  animals  recaptured  within  a  trapping  occasion  combined  across  occasions  (Residents) 
and  proportion  of  first-captures  only  (Transients)  for  three  species  of  small  mammals  at  Montezuma  Castle  Na- 
tional Monument, Yavapai  County,  Arizona.  *Log  likelihood  ratio  chi-square  value=4.837,  d.{.-l,p<0.05. 


19 


TECHNICAL  REPORT  NPS/NAUMOCA/NRTR-96/11 


Table  3.  Mean  (+SE)  percent  cover  of  vegetative  strata  and  ground  cover  by  habitat  type  (n=200  for 
each  habitat),  and  frequency  of  debris  piles  at  Montezuma  Castle  National  Monument,  Yavapai  County, 
Arizona.  Percent  cover  by  strata  were  analyzed  with  a  Wilcoxon  (Mann- Whitney)  signed  ranks  test. 
Debris  pile  frequency  differences  were  analyzed  using  a  G-test. 


Mean  percent  cover 


Vegetative  Strata 

Floodplain 

Trees 

28.02  ±  3.36 

Shrubs 

32.82  ±  2.30 

Forbs 

3.98  ±  0.87 

Perennial  grasses 

1.82  +  0.58 

Annual  grasses 

7.49  ±  1.28 

Debris  piles 

182 

Ground  Cover 

Bare  soil 

12.42  ±  1.41 

Gravel 

4.96  +  0.93 

Rock 

20.91  ±  2.25 

Litter 

48.75  +  2.38 

Mesquite 


43.35  ±  2.58  * 

29.62  ±  2.03 

0.03  ±  0.02  * 

6.46  ±  0.93* 

8.46  ±  1.24 

45** 


19.71  ±1.36* 

3.76  ±  0.68 

3.74  ±  0.26  * 

75.78  +  2.00  * 


*  p<0.05 

**  Log  likelihood  ratio  chi-square=l  10.7;  p<0.05 


total  of  three  species  of  perennial  grasses,  several 
species  of  annual  grasses,  two  species  of  herbs,  and 
1 2  different  shrub  species  were  encountered  in  mes- 
quite habitat. 

In  riparian  floodplain  habitat,  the  percentage 
of  perennial  grass  cover  was  an  important  environ- 
mental variable  as  indicated  by  the  length  of  its 
line  relative  to  other  lines  (Figure  9).  In  Figure  9, 
we  projected  a  perpendicular  dotted  line  from  P. 
eremicus  to  the  perennial  grass  axis  in  the  diagram 
and  the  endpoints  of  this  dotted  line  indicate  the 
relative  value  of  the  weighted  average  of  this  spe- 
cies with  respect  to  occurrence  of  perennial  grasses. 
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Therefore, we  inferred  that  P.  boylii  has  the  lowest 
weighted  average  with  respect  to  perennial  grass 
cover  (this  species  occurs  in  areas  with  high  litter 
content  and  annual  grass  cover),  N.  albigula  has  the 
second  lowest  value,  and  P.  eremicus  is  inferred  to 
have  the  highest  weighted  average.  Hence,  P. 
eremicus  was  captured  in  areas  with  a  high  percent- 
age of  perennial  grasses,  whereas  N.  albigula  and  P. 
boylii  were  captured  in  areas  with  lower  percent- 
ages of  perennial  grasses. 

In  floodplain  habitat,  we  found  percent  peren- 
nial grass  cover,  forb  cover,  tree  cover,  and  fre- 
quency of  trees  were  the  strongest  correlates  with 
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Table  4.    Average  percent  vegetation  cover  by  species  by  habitat  type  and  species  richness  by  habitat 
type.  Numbers  in  parentheses  are  standard  error  values.  (n=200  for  each  habitat  type). 


Mean  percent  cover 


Species 


Floodplain 

Perennial  grasses 

Bothriochloa  barbinodis 

0.09(0.05) 

Bouteloua  curtipendula 

1.60(0.57) 

Aristida  purpurea 

3.60(1.57) 

Muhlenbergia  porteri 

-0- 

Paspalutn  distichum 

0.14(0.08) 

Sitanion  hystrix 

-0- 

Sporobolus  airoides 

-0- 

Annual  grasses 

Bromus  spp. 

2.30(0.61) 

Unknown  annual  grass 

5.19(1.17) 

Herbs 

Castilleja  chromosa 

0.01(0.01) 

Datura  meteloides 

0.31(0.08) 

Eriodictyon  angustifolium 

0.14(0.11) 

Eriogonum  spp. 

-0- 

Linaria  dalmatica 

0.01(0.01) 

Lotus  mearnsii 

0.09(0.04) 

Solidago  wrightii 

0.35(0.11) 

Unknown  annual  herbs 

3.35(0.88) 

Shrubs 

Acacia  greggii 

0.71(0.34) 

Atriplex  canescens 

-0- 

Baccharis  sarothroides 

0.28(0.16) 

Berberis  haematocarpa 

0.67(0.32) 

Brickellia  califomica 

0.57(0.30) 

Ceratoides  lanata 

-0- 

Chihpsis  linearis 

28.66(2.18) 

Ephedra  viridus 

-0- 

Gutierrezia  sarothrae 

1.22(0.20) 

Larrea  divaricata 

-0- 

Lycium  pallidum 

-0- 

Marrubium  vulgare 

-0- 

Opuntia  leptocaulis 

0.05(0.05) 

Rhus  trilobata 

0.37(0.27) 

Yucca  elata 

-0- 

Ziziphus  obtusifolia 

0.01(0.01) 

Unknown  shrubs 

0.03(0.02) 

Trees 

Celtis  reticulata 

0.10(0.07) 

Fraxinus  pennsylvanica 

5.06(1.42) 

Juglans  major 

1.52(0.79) 

Juniperus  monosperma 

5.67(1.44) 

Platanus  wrightii 

9.86(1.96) 

Prosopis  velutina 

5.61(1.36) 

Total  cover 

77.6% 

Total  number  of  species 

29 

Mesquite 


-0- 

-0- 

-0- 
5.78(0.94) 

-0- 
0.65(0.18) 
0.03(0.03) 


8.37(1.24) 
0.09(0.09) 


-0- 

-0- 

-0- 
0.01(0.003) 

-0- 

-0- 

-0- 
0.02(0.02) 


10.26(1.47) 
3.18(0.60) 

-0- 
3.75(0.84) 

-0- 
0.07(0.05) 

-0- 
0.38(0.21) 
9.42(0.85) 
0.09(0.09) 
0.13(0.08) 
2.19(0.45) 

-0- 

-0- 
0.06(0.05) 
0.03(0.02) 
0.07(0.07) 


-0- 

-0- 

-0- 
0.38(0.37) 

-0- 
42.98(2.58) 

87.9% 
21 
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Figure  9.  The  distribution  of  three  species  of  small  mammals  in  the  riparian  floodplain  habitat  at  Montezuma 
Castle  National  Monument,  Yavapai  County,  Arizona.  Canonical  correspondence  analysis  (CCA)  ordination 
diagram  with  species  as  points  and  environmental  variables  as  lines  ending  in  points;  first  axis  is  horizontal, 
second  axis  is  vertical.  Shown  also  is  the  projection  of  the  small  mammal  point  labeled  PEER  onto  the  trajectory 
of  the  line  Perennial  grasses.  The  endpoint  of  this  projection  indicates  the  approximate  weighting  of  the  centers 
of  distributions  of  this  mammal  along  the  variable  "perennial  grasses",  Peromyscus  eremicus  being  found  in  habitats 
with  the  highest  percentage  cover  of  perennial  grasses.  The  small  mammal  species  are  (clockwise  from  left  top 
quadrant):  PEBO  =  Peromyscus  boylii,  PEER  =  P.  eremicus,  and  NEAL  =  Neotoma  albigula.  The  environmental 
variables  are:  Litter  =  percentage  substrate  covered  by  litter,  Annual  grasses  =  percentage  cover  of  annual  grasses, 
Gravel  =  percentage  substrate  covered  by  gravel,  Bare  soil  =  percentage  substrate  covered  by  bare  soil, 
Forbs  =  percentage  cover  of  herbaceouds  vegetation,  Shrubs  =  percentage  cover  of  shrub  species,  Shrub  #  = 
density  of  shrubs,  Rock  =  percentage  substrate  covered  by  rock,  Perennial  grasses  =  percentage  cover  of  perennial 
grasses,  Debris  =  frequency  of  debris  piles,  Trees  =  percentage  cover  by  trees,  Tree  #  =  density  of  trees,  and 
Sand  =  percentage  substrate  covered  by  sand. 


small  mammal  distributions.  Peromyscus  eremicus, 
displayed  in  the  lower  right  quadrant  of  the  dia- 
gram, was  captured  mainly  in  areas  with  high  per- 
centage of  cover  by  perennial  grasses  and  where 
substrate  was  fairly  rocky.  Peromyscus  boylii,  in  the 
left  upper  quadrant,  was  captured  mainly  in  areas 
with  high  annual  grass  cover,  high  fitter  contents, 
high  percent  cover  by  forbs,  and  where  substrate 
was  mostly  exposed  soil.  Neotoma  albigula,  displayed 
in  the  lower  left  quadrant,  occurred  in  areas  with 
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high  fitter  contents,  and  high  tree  frequency.  For 
this  habitat  type,  100.0%  of  the  cumulative  vari- 
ance of  the  species-environment  relationship  was 
accounted  for  by  the  first  two  axes  (Table  5). 

Microhabitat  separation  among  these  three  spe- 
cies in  the  floodplain  was  significant  (Table  5)  .The 
99  random  data  sets  generated  by  the  Monte  Carlo 
permutation  test  of  trap  stations  yielded  a  lower 
eigenvalue  with  a  p- value  less  than  0.05.  There- 
fore, differences  in  vegetative  strata  and  habitat 
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Table  5.  Summary  of  canonical  correspondence  analysis  (CCA)  and  Monte  Carlo  permutation  test  in 
desert  riparian  floodplain  habitat  at  Montezuma  Castle  National  Monument, Yavapai  County,  Arizona. 
Results  are  from  a  CCA  of  three  small  mammal  species  and  1 3  environmental  variables  shown  in  the 
ordination  diagram  of  Figure  9. 


Axes 

1 

2 

3 

4 

Eigenvalues: 

.130 

.064 

.521 

.351 

Species-environment  correlations: 

.518 

.330 

.000 

.000 

Cumulative  percentage  variance 

of  species  data: 

12.2 

18.1 

67.0 

100.0 

of  species-environment 

relation: 

67.1 

100.0 

.000 

.000 

Sum  of  all  canonical  eigenvalues: 

Monte  Carlo  permutation  test  for  significance  of  first  canonical  axis  (1): 
Eigenvalue  =  0.13 
F-ratio  =  14.68 
*  P-value  =  0.04 


1.066 


structure  found  at  a  trap  station  significantly  de- 
termine what  species  will  occur  there. 

In  mesquite  habitat,we  found  P.  eremicus  also 
occurred  in  areas  with  high  percentage  of  peren- 
nial grasses,  and  areas  dense  with  shrubs  and  her- 
baceous species  (Figure  10).  This  species  was  nega- 
tively associated  with  tree  cover,  tree  frequency, 
litter  content,  and  debris  piles.  As  in  floodplain 
habitat,  N.  albigula  occurred  most  often  with  debris 
piles,  high  percent  tree  cover,  tree  density,  and  high 
percent  litter  content.  Perotnyscus  boylii,  in  the  up- 
per right  quadrant,  occurred  in  areas  with  high 
percentage  of  annual  grasses,  and  where  the  sub- 
strate is  mosdy  rocky.  For  this  habitat,  100%  of  the 
cumulative  variance  of  the  species-environment 
relation  can  be  explained  by  the  first  two  axes  (Table  6). 

Microhabitat  separation  among  these  three  spe- 
cies in  the  mesquite  habitat  was  also  significant 
(Table  6).  The  Monte  Carlo  permutation  test  of 
trap  station  again  yielded  a  lower  eigenvalue  for 
the  first  canonical  axis  indicating  significant  mi- 
crohabitat segregation  among  species  within  a  habi- 
tat type. 

Across  habitat  types,  the  difference  between 
microhabitats  where  P.  boylii,  P.  eremicus,  and  N. 
albigula  were  caught  and  those  where  these  species 


were  never  caught  was  significant  (F—  5.98,Wilks' 
Lambda=0.769,  d.f.  =  10,  199,  ;><0.05;  F=2.09, 
Wilks*  Lambda-0.899,  d.f.=10, 185./K0.05;  and, 
F=3.55,  Wilks'  Lambda=0.761,  d.f.  =  10,  113, 
p<0.05  respectively)  (Appendix  B).  Results  from 
the  univariate  F-tests  showed  that  Perotnyscus  boylii 
was  captured  at  trap  stations  where  the  percent 
shrub  cover  was  significantly  higher  (F=  11.14, 
d.f.=l,/><0.05),  where  frequency  of  debris  piles 
was  higher  (G=66.34,  d.f.=l,p<0.05),  and  where 
percent  rocky  (F=26.49,d.f.=l,p<0.05)  and  sandy 
(F=13.18,  d.{.=l,p<0.05)  substrates  were  signifi- 
candy  higher  than  randomly  available,  but  not  used, 
trap  stations.  Perotnyscus  boylii  was  also  captured  at 
stations  with  less  annual  grass  cover  (F=5.91,  d.f.=l, 
p<0.05),  less  percent  cover  by  litter  (F=20.82, 
d.f.  =  l,  p<0.05),  and  lower  frequency  of  trees 
(G=9.86,  di>l,jK0.05). 

Perotnyscus  eremicus  was  captured  at  trap  stations 
with  greater  percent  shrub  cover  (F=8.64,  d.f.=l, 
p<0.05),  a  greater  frequency  of  shrubs  (G=3.80, 
d.f.=l,  p<0.05),  significantly  rockier  substrate 
(F=4.91,  d.f.=l,  p<0.05),  and  less  exposed  soil 
(F=4.92,  d.f.=l,/K0.05). 

Neotoma  albigula  was  captured  in  areas  -with 
greater  tree  cover  (F=6.38,  d.f.=l,p<0.05),  higher 
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Figure  10. The  distribution  of  three  species  of  small  mammals  in  mesquite  habitat  at  Montezuma  Castle  National 
Monument, Yavapai  County,  Arizona.  Canonical  correspondence  analysis  (CCA)  ordination  diagramwith  spe- 
cies as  points  and  environmental  variables  as  lines  ending  in  points;  first  axis  is  horizontal,  second  axis  is  vertical. 
The  small  mammal  species  are  (clockwise  from  left  top  quadrant):  PEER=  Peromyscus  eremicus,  PEBO=  P. 
boylii,  and  NEAL=  N.  albifula.  The  environmental  variables  are:  Perennial  grasses=  percentage  cover  by  perennial 
grasses,  Gravel=  percentage  substrate  covered  by  gravel,  Forbs=  percentage  cover  by  herbaceous  vegetation,  Bare 
soil=  percentage  substrate  covered  by  bare  soil,  Annnual  grasses=  percentage  cover  by  annual  grasses,  Rock= 
percentage  substrate  covered  by  rock,  Litter=  percentage  substrate  covered  by  Utter,  Debris  =  frequency  of  debris 
piles,  Trees=  percentage  cover  by  trees,  Tree  #=  density  of  trees,  Shrub  #=  density  of  shrub  species,  Shrubs= 
percentage  cover  of  shrub  species. 


frequency  of  debris  piles  (G=28.81,  d.f.  =  l, 
p<0.05),  and  higher  frequency  of  trees  (G=7.29, 
d.f.=l,jK0.05).This  species  was  also  captured  in 
areas  with  less  exposed  soil  (F=  15.52,  d.f.=l, 
^<0.05)  and  less  cover  by  perennial  grasses  (F=8.66, 
d.f.=l,/K0.05). 

Within  the  floodplain  habitat,  there  were  no 
significant  differences  in  percent  cover  of  vegeta- 
tive strata  for  traps  used  and  traps  randomly  avail- 
able within  a  habitat  type  (Appendix  C).  Structural 
characteristics  of  the  habitat  did  differ  between  used 


and  randomly  available  trap  stations  for  P.  boylii  and 
N.  albigula  with  P.  boylii  associated  with  a  higher 
frequency  of  debris  piles  (G=17.81,  d.f.  =  l, 
p<0.05),  and  N.  albigula  also  associated  with  more 
debris  (G=17.28,d.f.=l,/K0.05)  and  higher  den- 
sity of  trees  (G=  13.70,  d.f.=l,p<0.05). 

Within  the  mesquite  habitat,  there  were  also 
no  significant  differences  in  vegetative  strata  be- 
tween used  and  unused  traps  (Appendix  D). Again, 
N.  albigula  was  associated  with  more  debris  (G=4.97, 
d.f.=l,/j<0.05),  a  structural  characteristic. 
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Table  6.  Summary  of  canonical  correspondence  analysis  (CCA)  and  Monte  Carlo  permutation  test 
results  in  mesquite  habitat  at  Montezuma  Castle  National  Monument,Yavapai  County,  Arizona.  Results 
are  from  a  CCA  of  three  small  mammal  species  and  12  environmental  variables  shown  in  the  ordination 
diagram  of  Figure  10. 


Axes 


4 


Eigenvalues: 

Species-environment  correlations: 
Cumulative  percentage  variance 

of  species  data: 

of  species-environment 

relation: 


.236 
.594 

18.4 

62.0 


.145 
.487 

29.8 

100.0 


.475 
.000 

66.9 

.000 


.422 
.000 

100.0 

.000 


Sum  of  all  canonical  eigenvalues: 


1.277 


Monte  Carlo  permutation  test  for  significance  of  first  canonical  axis  (1): 
Eigenvalue  =  0.24 
F-ratio  =  12.67 
*  P-value  =  0.02 
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CHAPTER  4 
DISCUSSION 


Habitat  quality  differences  in  a  heterogeneous 
landscape  can  result  in  a  source-sink  population 
structure  with  net  emigration  of  individuals  from 
better  habitats  (sources)  combined  with  net  immi- 
gration into  poorer  habitats  (sinks)  (Pulliam  1988; 
Pulliam  and  Danielson  1991).  This  source-sink 
population  structure  manifests  differences  in  per- 
formance and  reproductive  success  of  individuals 
encountering  the  source  and  sink  habitats  (Kawecki 
1995).  Hence,we  define  high  quality  habitats,  or 
sources,  to  be  those  areas  that  afford  conditions 
necessary  for  relatively  successful  survival  and  re- 
production over  long  periods  when  compared  with 
other  similar  environments.  Conversely,  marginal 
habitats,  or  sinks,  support  individuals,  but  their  rates 
of  survival  and  reproduction  are  low  relative  to  high 
quality  habitats.  Marginal  habitats  are  usually  suit- 
able for  occupancy  for  short,  or  intermittent,  pe- 
riods (Morrison  et  al.  1992).  In  a  purely  social  sense, 
dispersal  sinks  can  also  develop  if  social  interac- 
tions prevent  subordinate  individuals  from  enter- 
ing into,  or  remaining  in,  high  quality  habitats  (Van 
Home  1982, 1983).  So,  sinks  can  also  be  thought 
of  as  marginal  areas  to  which  surplus  individuals 
disperse  due  to  competitive  exclusion  from  opti- 
mal habitats  (Lidicker  1976). 

Source  and  sink  dynamics  within  a  population 
have  been  documented  in  a  variety  of  taxa,  e.g. 
Bonasa  bonasia  (hazel  grouse)  (Beshkarev  et  al. 
1994),  Peromyscus  maniculatus  (deer  mouse)  (Van 
Home  1982),  Tamias  townsendii  (Townsend's  chip- 
munk), Glaucomys  sabrinus  (northern  flying  squir- 
rel), Zapus  trinotatus  (Pacific  jumping  mouse)  (Doyle 
1990),  sphingid  moths  (Janzen  1986),  Rangifer 
tarandus  (caribou)  (Bergerud  1988),  and  Haleaetus 
leucocephalus  (bald  eagles)  (Swenson  et  al.  1986). 
These  studies  show  that  the  source-sink  phenom- 
enon is  widespread,  but  difficult  to  determine  con- 
clusively within  a  study  site  unless  an  excess  of  local 
reproduction  over  immigration  can  be  shown,  and 
active  juvenile  dispersal  is  observed  (Andersen 
1994). 

Since  source  habitats,  by  definition,  are  net 
exporters  of  individuals,  and  are  presumed  of  higher 


quality  than  adjacent  habitats,  we  made  several  pre- 
dictions about  the  small  mammal  community  struc- 
ture within  a  source.  These  included  a  greater 
population  density,  greater  species  diversity,  larger- 
sized  adults,  and  a  greater  proportion  of  the  popu- 
lation in  reproductive  condition  within  a  source 
for  a  given  species.  In  addition,  survival  of  indi- 
viduals should  increase  in  these  high  quality  habi- 
tats, number  of  transient  animals  should  be 
minimized,  and  most  importantly,  there  should  be 
evidence  of  dispersal  from  source  to  sink. 

Riparian  floodplains  are  generally  considered 
more  productive  than  adjacent  ecosystems  because 
of  their  unique  hydrologic  conditions  (Mitsch  and 
Gosselink  1993),  and  this  is  particularly  true  in  the 
arid  southwest  (Johnson  1979).  Periodic  flooding 
contributes  to  a  higher  productivity  by  supplying 
an  adequate  water  supply  for  the  vegetation,  sup- 
plying nutrients  and  altering  soil  chemistry  (in- 
creasing nitrification,  sulfate  reduction,  and  nutrient 
mineralization),  and  creating  a  more  oxygenated 
root  zone  while  also  "flushing"  waste  products. 

For  these  reasons,  we  predicted  that  on  a  re- 
gional scale  the  floodplain  would  be  the  source  for 
several  small  mammal  species.  But  the  data  that  we 
collected  for  floodplain  and  mesquite  habitats  at 
Montezuma  Castle  did  not  support  this  prediction. 
We  also  predicted  a  more  diverse  small  mammal 
community  in  the  floodplains  since  a  habitat  of 
high  quality  should  be  able  to  support  a  more  di- 
verse small  mammal  assemblage.  But  we  found  spe- 
cies diversity  greater  in  the  mesquite,  due  to  the 
presence  of  burrowing  small  mammals  and  habi- 
tat-specific species.  For  instance,  Dipodomys  ordii 
(Ord's  kangaroo  rat)  was  captured  only  in  mes- 
quite where  the  substrate  included  a  greater  amount 
of  sandy,  gravelly,  and  friable  soil.  The  floodplain 
areas  also  contained  friable  soil,  but  the  occurrence 
of  this  substrate  was  patchy  due  to  fragmentation 
by  stream  channels.  In  addition,  Reithrodontomys 
megalotis  (western  harvest  mice)  was  captured  only 
in  mesquite  habitats. 

Although  abundance  of  P  boylii  was  greater  in 
the  floodplain,  male  adults  were  not  significandy 
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larger,  nor  was  there  an  increase  in  the  number  of 
juveniles  and  subadults  in  comparison  to  mesquite 
habitat.  In  addition,  reproductive  activity  was  not 
greater  in  the  floodplain,  although  non-significance 
for  these  latter  two  criteria  may  be  a  reflection  of 
low  sample  sizes  in  mesquite  habitat.  Finally,  P.  boylii 
did  not  appear  to  survive  significantly  longer  in 
floodplain  over  mesquite  habitats. 

Peromyscus  boylii  was  more  transient  in  mesquite. 
Animals  in  poor  quality  habitat  or  in  a  sink  habitat 
might  continue  searching  for  better  habitat  and 
would  gain  little  advantage  from  the  long-term 
maintenance  of  a  stable  home  range,  whereas  sta- 
bility might  well  be  important  for  both  dominants 
and  subordinates  in  high  quality  habitat  (Van  Home 
1981). 

Peromyscus  eremicus  population  size  also  tended 
to  be  larger  in  the  floodplain,  but  again  results  from 
the  other  criteria  were  not  conclusive.  There  was 
some  indication  that  mesquite  was  a  source  habitat 
for  N.  albigula.  Neotoma  albigula  population  size  was 
generally  larger  in  the  floodplain,  but  males  were 
significantly  larger  in  the  mesquite,  and  there  were 
more  juveniles  and  subadults  in  the  floodplain 
which  would  indicate  mesquite  is  a  higher  quality 
habitat  for  this  species.  Sinks  can  exhibit  larger 
population  sizes  relative  to  sources  (Pulliam  1988; 
Van  Home  1983). 

We  could  not  conclude  from  our  data  that  ri- 
parian floodplain  habitat  along  Wet  Beaver  Creek 
acted  as  a  source  habitat  and  that  mesquite  acted  as 
a  dispersal  sink  for  any  small  mammal  species. The 
fact  that  our  data  do  not  support  source-sink  dy- 
namics may  suggest  three  things.  First,  the  spatial 
scale  and  the  number  of  sampling  replications  of 
the  study  site  were  both  relatively  small.  On  a  land- 
scape level,  riparian  floodplain  may  act  as  a  source 
for  P.  boylii,  but  within  our  study  site,  we  were  not 
able  to  detect  this  pattern.  We  also  did  not  sample 
true  "upland"  habitat.  Second,  since  population 
densities  appeared  to  be  low  during  the  study,  in- 
terspecific competition  or  interference  among  in- 
dividuals may  not  have  been  sufficiently  strong  to 
force  emigration  into  inferior  habitats.  Third,  the 
temporal  scale  of  this  study  may  not  have  been 
sufficient  to  detect  a  source  and  sink  habitat. 


Although  we  were  not  able  to  distinguish  a 
source  from  a  sink  for  any  species,  the  data  that  we 
collected  on  habitat  use  provides  important  infor- 
mation on  habitat  selection  and  habitat  separation 
among  the  species.  The  distribution  of  individuals 
among  habitats  may  be  largely  determined  by  habi- 
tat selection  (Pulliam  and  Danielson  1991).  At  low 
population  densities,  habitat  selectivity  is  enhanced, 
but  at  high  densities,  habitat  selectivity  diminishes 
(Rosenzweig  1991).  Habitat  selectivity  by  P.  boylii, 
P.  eremicus,  and  N.  albigula  was  significant  between 
floodplain  and  mesquite  habitats  indicating  that 
selection  of  habitat  was  probably  the  most  impor- 
tant mechanism  structuring  the  small  mammal 
communities  within  Montezuma  Castle  NM. 

Peromyscus  boylii  is  a  widely  distributed  species 
living  in  a  variety  of  habitats  in  Arizona.  This  spe- 
cies is  a  large-sized  Peromyscus  and  frequently  climbs 
in  and  through  trees  with  ease  (Hoffmeister  1986). 
During  this  study,  we  found  P.  boylii  to  be  more 
abundant  in  the  floodplain  where  it  selected  mi- 
crohabitats  with  more  debris,  sandy  soils,  and  a 
rocky  substrate.  In  the  Huachuca  Mountains  of 
Arizona,  Hoffmeister  and  Goodpaster  (1954)  also 
found  P.  boylii  to  be  the  most  common  mammal  in 
tree  covered  areas  with  rocky  substrate  and  heavy 
undergrowth  and  in  riparian  or  wash  habitat  along 
streams.  Goodwin  and  Hungerford  (1979)  also 
found  P.  boylii  in  high  densities  along  rocky  slopes 
and  brush-rock  slopes.  Association  of  P.  boylii  with 
floodplain  habitats  probably  also  relates  to  the  her- 
baceous vegetation  in  these  habitat  types.This  spe- 
cies has  a  flexible  diet,  but  the  main  component  is 
herbaceous  vegetation  (Goodwin  and  Hungerford 
1 979) .  Percent  cover  by  herbaceous  vegetation  was 
significantly  higher  in  the  floodplain  habitat  than 
in  the  mesquite  in  our  study,  which  could  partially 
explain  why  this  species  is  more  abundant  in  the 
floodplain. 

Peromyscus  eremicus  is  the  most  typical  and  wide- 
spread desert-dwelling  member  of  the  genus 
Peromyscus  (King  1968)  .This  species  is  found  among 
and  around  cacti,  creosote,  woodpiles,  rocks  and 
rocky  slopes,  chaparral,  and  sandy  flats,  and  they 
often  reside  in  the  abandoned  burrows  of  other 
mammals  or  in  parts  of  woodrat  piles  (Hoffmeister 
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1 986) .  Peromyscus  eremicus  at  Montezuma  Castle  also 
tended  to  be  more  abundant  in  the  floodplain.This 
species  selected  microhabitats  with  significantly 
more  shrub  cover,  significantly  higher  density  of 
shrubs,  and  a  rockier  substrate  than  randomly  avail- 
able microhabitats.  In  coastal  sage  scrub  habitat  of 
California,  frequency  of  capture  of  P.  eremicus  was 
also  greater  in  areas  of  dense  vegetation  and  shrubs 
(Veal  and  Caire  1979). 

Neotoma  albigula  is  common  in  a  variety  of  habi- 
tats throughout  Arizona,  building  large  nests  in  and 
around  cacti,  shrubs  such  as  mesquite  or  acacia,  or 
rock  piles.  In  mesquite  and  catclaw  vegetation  as- 
sociations, this  species  feeds  on  the  bark  of  these 
shrubs  frequendy  climbing  into  them  to  debark 
certain  branches  (Hoffmeister  1986)  .We  found  that 
N.  albigula  selected  microhabitats  with  greater  tree 
cover,  increased  tree  density,  and  greater  amount 
of  debris  than  microhabitats  that  were  available, 
which  follows  habitat  utilization  patterns  reported 
in  the  literature.  Macedo  and  Mares  (1988)  also 
found  N.  albigula  in  areas  with  large  quantities  of 
debris  and  dead  trees.  Large  quantities  of  debris 
provide  sufficient  cover  for  house  construction. 
Olsen  (1973)  hypothesized  that  shelter  site  selec- 
tion by  this  species  is  based  on  the  quantity  of 
ground  level  vegetation  and  debris  available  for 
cover. 

Not  only  were  the  above  species  of  small  mam- 
mals selecting  certain  habitat  characteristics 


nonrandomly,  but  we  also  found  evidence  of  sig- 
nificant microhabitat  separation  among  these  spe- 
cies. In  floodplain  habitat,  N.  albigula  was  associated 
with  more  trees  and  tree  cover,  P.  eremicus  with 
rockier  substrates  and  increased  cover  by  perennial 
grasses,  while  P.  boylii  was  more  highly  associated 
with  herbaceous  vegetation,  exposed  soil  substrate, 
and  cover  by  annual  grasses.  In  mesquite,  P  boylii 
was  associated  with  more  annual  grass  cover  and  a 
rockier  substrate,  P.  eremicus  with  perennial  grass 
cover,  herbaceous  vegetation,  and  gravelly  soils, 
while  N.  albigula  was  again  associated  with  more 
trees  and  tree  cover.  Neotoma  albigula  was  also  asso- 
ciated more  significantly  with  debris  and  litter.  Co- 
existence of  small  mammal  species  within  the  same 
area  is  facilitated  by  microhabitat  separation  (Doyle 
1987;  Dueser  and  Shugart  1978),  and  microhabi- 
tat separation  has  been  shown  to  be  a  prominent 
feature  in  community  structure  in  several  studies 
in  the  southwest  (Holbrook  1979;Rosenzweigand 
Winakur  1969;  Thompson  1982). 

Both  the  importance  of  habitat  selection  and 
the  role  that  microhabitat  plays  in  structuring  com- 
munities is  -well  documented  (Adler  1988;  Price 
1978).  Small  mammal  community  structures  in 
riparian  floodplain  and  mesquite  habitats  at 
Montezuma  Castle  National  Monument  in  cen- 
tral Arizona  may  depend  on  several  factors,  but 
habitat  selection  and  microhabitat  separation  ap- 
pear to  be  the  most  influential  mechanisms. 
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As  the  nation's  principal  conservation  agency,  the  Department  of  the  Interior  has  responsibility  for  most  of 
our  nationally  owned  public  lands  and  natural  and  cultural  resources.  This  includes  fostering  of  our  land  and 
water  resources,  protecting  our  fish  and  wildlife,  preserving  the  environmental  and  cultural  values  of  our 
national  parks  and  historical  places,  and  providing  for  enjoyment  of  life  through  outdoor  recreation.  The 
department  assesses  our  energy  and  mineral  resources  and  works  to  ensure  that  their  development  is  in  the 
best  interests  of  all  our  people.  The  department  also  promotes  the  goals  of  the  Take  Pride  in  America  campaign 
by  encouraging  stewardship  and  citizen  responsibility  for  the  public  lands  and  promoting  citizen  participation 
in  their  care.  The  department  also  has  a  major  responsibility  for  American  Indian  reservation  communities 
and  for  people  who  live  in  island  territories  under  U.S.  administration. 
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